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FAR ULTRAVIOLET EXCITATION PROCESSES IN COMETS 
P. D. Feldman, C. B. Opal, R. R. Meier and K. R. Nicolas 


Introduction 


The recent observations of atomic oxygen and carbon 
in the far ultraviolet spectrum of Comet Kohoutek (1973f) 
(Feldman e^ al. 1974; Opal e^ al. 1974) have demonstrated 
the existence of these atomic species in the cometary 
coma. However, in order to identify the source of their 
origin, it is necessary to relate the observed ultraviolet 
flux to the atomic production rate. Assuming the only 
excitation mechanisms allowed are those produced by 
resonance scattering and fluorescence of solar ultra- 
violet radiation, the problem reduces to finding the 
emission rate factors (g-factors) as a function of the 
heliocentric comet velocity. Since the widths of the 
solar emission lines are smaller than the maximum helio- 
centric Doppler shift, given by 

M G 

r = (-2-) = 21.06 km sec"^ , 

max q 

where q is the perihelion distance in A.U., it is 
necessary to consider the detailed multiplet structure 
of the transition, the solar line shape and the relaxation 
of excited fine structure levels. 

Analyses of the observed 01 A1304 and Cl M657 A 
multiplets have been carried out using high resolution 
solar spectra obtained from the ATM solar spectrograph 
aboard Sky lab. In addition, we have examined the 
possibility of observing ultraviolet fluorescence from 
molecules such as CO (which may be the parent molecule 
of atomic carbon) and H 2 , as well as resonance scattering 
either from atomic ions'^for which there are strong 
corresponding solar lines (CII) or from atoms for which 
there is an accidentical wavelength coincidence (SI) . 

The scattering of solar Lycx from atomic hydrogen has been 
discussed in detail by Keller (1973) and Meier (1974) and 
will not be considered here. 


Emission Rate Factor 


The emission rate factor, which is the probability 
that an atom or molecule will resonantly scatter a solar 
photon of wavelength A into 4ir sr in unit time, is given 
by Barth ( 1969 ) as: 

2 

^A^ photons sec ^ mol~^ 
me 


where f^ is the absorption oscillator strength and ttF^ 
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the solar flux per unit wavelength. For resonance 
fluorescence in molecules, the relative transition 
probabilities for downward transitions must also be taken 
into account. The emission rate per unit volume is then 
simply related to the density of the scattering species 
by 

- 1 -3 

^i Photons sec cm 


In the case of cometary emission, the large helio- 
centric velocity of a typical comet requires the use 
of the solar flux ttF^i , where 

I r V 

>> = ^(1 - 1 ) , 


so that g^ = g^(r) = g^(r, q) . 


Solar Ultraviolet Fluxes 


The Naval Research Laboratory solar spectrograph 
(S082B) on the Apollo Telescope Mount (ATM) of Skylab 
photographed the solar spectrum from 970 A - 4000 A. It 
provided the necessary high spectral resolution (0.06 A) 
to yield data necessary for the cometarv interpretation. 
The spectrograph slit averaged over a 2" by 60” area on 
the sun. In order to estimate the solar intensity from 
the total disk it was necessary to include the effects 
from limb brightening and from the increased contribution 
of active regions which introduces a variable intensity 
component. Spectra of these phenomena were available 
from the large number of observations taken within the 
Joint Observing Program. 

Densitometry of the Kodak 104 film gave the density 
versus position information which was used to determine 
the relative intensity versus wavelength of the solar 
spectrum. "Spectra with 40, 160, and 640 sec exposure 
times were used to construct a relative H-D curve for 
each wavelength region of interest. The conversion from 
density to relative intensity for each of the plates in 
the exposure sequences yielded a maximum variation in 
relative line shape of about i’15%, the random error being 
mostly due to grain clumping. The wavelength determina- 
tion was made by using from 5 to 9 standard lines 
(Sandlin, 1974) within each wavelength interval. The 
estimated absolute error is i’0.02 A which translates into 
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i’lOTo intensity error at the steepest part of the line 
profile. 

The absolute intensity was determined independently 
of the ATM calibration in two ways.* The first method 
was to match the continuum intensity level at 1650 A 
with a value determined from a 1971 rocket flight 
(Brueckner and Nicolas, 1972). This gave an absolute 
intensity scale for the Cl emission feature at 1656 A. 

The ms error is about ±307o. At shorter wavelengths the 
absolute scale was determined by comparing the total 
relative line intensities with those of OSO VI (Dupree 
et al. 1973). The error for this method is approximately 
+r0D^ -50%. 

The total flux averaged over the solar surface and 
its time variation was estimated from spectra of quiet 
solar limb scans and of active regions . The intensity 
of the emission lines is fairly constant- (neutral) over 
the solar disk. The maximum error introduced by assuming 
that they have neutral limb brightening is less than i’10%. 
The active region enhancement over the quiet region 
intensity is approximately a factor of 2 to 3 for Cl and 
a factor of 10 to 20 for 01 and CII. Since the area of 
the disk covered by active regions can be as high as 107o 
(De Jager 1961), the total solar flux during the solar 
cycle could vary by 307o for Cl and up to a factor of 2 
to 3 for 01 or CII. Short tern fluctuations (on the 
order of two weeks) within a solar cycle could also be 
of comparable magnitude. 

Cl M657 A 


The carbon multiplet at 1657 A provides an 
excellent example of the dependence of g^ on radial 
velocity. Carbon, long known as a major constituent 
of cometary molecules and radicals, was not detected in 
atomic form until the two recent rocket ultraviolet 
observations of Comet Kohoutek. There are six components 
of the multiplet, as shown in Fig. 1, resulting from the 
fine structure triplet nature of both the excited and 
ground state. An additional consideration in the evalu- 
ation of the g-factors in cometary emission is the question 
of relative population of ground state fine structure 
levels. We assume that they are initially populated 
according to a Boltzmann distribution characteristic of 
some temperature Tq of the parent molecule. However, 
since the coma becomes non-collisional very close to 
the nucleus, radiative transitions J=2^ J=1 and 
J=1 J=0 will tend to cool the atom in a time of the 


The inflight ATM calibration is not yet finished, but 
will be completed and published shortly. 
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Fig. 1, High resolution spectrum (aX = 0. 06 A) of the sun in the region of the Cl multiplet at 1657 A, observed 
300" in from the solar limb. 



order T = ^ii” > '^here is the Einstein coefficient 
for the transition. The efrect of cooling is significant 
only if the A^j, values are greater than the photoionization 
rate which depends on heliocentric distance as r"2. 

For the carbon ground state, A^q = 7.9 x 10~8 sec"l- and 
At]^ = 2.7 X 10"7 sec“l-, which are both less than the value 
of J4 at 1 A.U. of 4.0 X 10“6 sec"l (Feldman et al. 1974). 
The distance at which the cooling time TT^ and”Tonization 
lifetime are equal, r^,, is given in Table I for the 
most important species. For cases in which the 

variation of the relative populations in time will be 
determined primarily by the effect of optical pumping 
produced by the resonance scattering of solar radiation, 
which as we have noted will also vary with both r and r. 

In most cases of interest, the resonance scattering time, 

~ is less than Xf., so that the optical pumping 

effect is not negligible. 

To first order we regard the temperature in the 
coma as constant and show in Figs. 2 and 3 the variation 
of g^ with r and r. All g-factors are calculated for 
r = 1 A.U. The values of r corresponding to the observa* 
tions of Comet Kohoutek are indicated in For 

planetary atmospheres the value of g for r = 0 is to be 
used. 


01 A 1304 A 


The oxygen triplet in the solar spectrum has 
recently been observed with sufficient resolution to 
permit the evaluation of the absorption at the line 
center due to terrestrial oxygen in the upper atmosphere. 
An example of one of the lines (after correcpion for 
atmospheric absorption) is shown in Fig. 4. A Doppler 
shift of +0.17 A corresponding to a velocity of +40 km 
sec”^ is sufficient to completely shift the cometary 
absorption wavelength completely off of the solar line. 
Thus the observed oxygen emission (Feldman £t al . 1974; 
Opal £t al. 1974) must be due to either resonance 
scattering from the solar continuum, which is very weak, 
or another mechanism. Fig. 5 illustrates a Bowen 
fluorescence mechanism in which the state is populated 
via the 3[)_3p transition at 1025.77 A which is nearly 
degenerate with the solar Ly P line of HI at 1025.72 A. 

The g-factors for the entire raultiplet for the two 
processes are illustrated in Fig. 4. Resonance scattering 
from the solar continuum gives a g-factor of 4 x 10”8 
sec"^ atom'l, independent of wavelength. For the case of 
oxygen, cooling via the fine structure transitions at 
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Tab le I 


Species 

01 

Cl 


COOLING AND IONIZATION TIMES 


T^(sec) 

T. (sec at 1 A.U.) 

r (A,U J 

10^ 

4.0 X 10^ 

— c 

0.05 

107 

2.5 X 10^ 

6.3 

4 X 10^ 


V «• 


CII 
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Fig. 2. Dependence of the Cl 1657 A g-factor on helio- 
centric velocity. 
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Dependence of the Cl 1657 A g-factor on hello 
centric distance for Comet Kohoutek (1973f). 
The arrows indicate direction of motion. 
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ig. 4. High resolution profile of the solar 01 1302.2 A 
line. The 01 g- factor is shown for both 
resonance scattering and fluorescence excited 
by solar LyP as a function of heliocentric 
velocity. The arrows correspond to the 
observations of Feldman et al . (1974) and 
Opal et al. (1974). 
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63m. and 147M- is quite rapid with 't However, 

for any reasonable range of equivalent temperatures, 
since the individual lines do not overlap, there is 
little variation in the total multiplet intensity. For 
the Bowen mechanism we assume that all of the atoms 
are in the lowest (J=2) level. 

While the intensity due to Bowen fluorescence is 
more than an order of magnitude smaller than that due 
to direct resonance scattering, for comet velocities 
^40 km sec"i it is the dominant source of cometary 
A 1304 emission. Using the Bowen g- factor, the oxygen 
production rate derived by Feldman e_t al. (1974) for 
Comet Kohoutek on 5.1 January 1974 was found to satisfy 
the requirement Qq ~ Qqh ~ ^HoO expected if 

photodissociation of water is the dominant source of 
atomic oxygen. At smaller heliocentric velocities the 
comet might be expected to be extremely bright at 1304 A 
and in fact could provide, near perihelion, a measurement 
of the solar 1304 A multiplet profile free of the effects 
of terrestrial oxygen absorption since the cometary 
absorption width is small compared with the width of 
the solar emission profile. The role of photoelectron 
excitation has been neglected since it is assumed that 
the coma becomes collisionless at radii greater than 
10^ km in which case the 1304 A image would appear to 
be less than 1 minute of arc contrary to the 10 arc min 
image reported by Opal ^ al . (1974). 

Because of the steep slope in the solar line shape, 
corresponding to Doppler velocities of the order of i"25 
km sec“i, an appreciable intensity asymmetry (Greens tein, 
1958) can result. To illustrate the effect of a rapidly 
varying solar line profile, we have computed atomic 
oxygen 1304 A intensities from Comet Kohoutek under some- 
what idealized conditions. For resonance scattering, 
one must compute the probability that a photon of a given 
frequency and direction (from the sun) will be scattered 
at a new frequency in the direction of the observer. The 
procedure for doing this is described by Meier (1974). 

At each point in the atmosphere this probability function 
must be integrated over the solar line, the velocity 
distribution of oxygen, and the emitted line profile. 

The resulting volume emission rate must then be integrated 
along the line of sight to compute the column emission 
rate. Assuming a radial outflow density distribution 
and a radial maxwell- Boltzmann velocity distribution 
(Keller, 1973; Bertaux £t al . , 1973), we have computed 
the expected Greens tein effect . The idealized observation 
for the comet is taken to be at a 90* sun-comet-earth 
angle. Velocities of 0 and -22.5 km sec'l relative to 
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the sun are used. The mean velocity of oxygen atoms is 
taken to be 1 and 2 km sec“i both with a production 
rate of 10^° atoms sec“l sr"^ and a g-factor appropriate 
to 1 A.U. The solar line shape in Fig. 4 was used. The 
line-center point of the comet rest frame is at r=0 in 
Fig. 4. As discussed above, all atoms are assumed to e 
be in the ground state so that although only the 1302.2 A 
line participates in excitation, all three lines will be 
present in emission. 

The results of this calculation are shown in Fig. 6. 
The larger asymmetry found for the 2 km sec"l outflow 
velocity is due to the broader cometary absorption 
coefficient ^overlapping a larger portion of the solar 
line. The r = -22^5 km sec“l isophotes have both been 
normalized to the r = 0 case (a factor of 1.275) for 
purposes of comparison in the figure. Thus it is clear 
that the degree of asymmetry upsun and downsun can yield 
information about the mean velocity of the scattering 
gas when the solar line varies significantly over the 
absorption line. 


CII A1335 A 

An anti-solar tail at least 5 x 10 km long and 
about 5 X 105 km wide was observed in the 1250-1800 A 
band before perihelion at r = .182 A.U. with the electro- 
graphic camera on Sky lab 4 (Page, 1974). This feature 
should also have been observable by the rocket instru- 
ments, assuming an r"^ power law dependence of brightness 
on sun-comet distance, however the rocket images in that 
band^of Opal et al. (1974) showed only a circular coma, 
attributed maTnly to 1304 A oxygen emission, and the 
spectrometer of Feldman ^ al. (1974) observed down-sun 
of the comet for over 30 sec without detecting any 
emissions from the comet. 

The tail observed from Skylab could be attributed 
to dust, a neutral species, or an ion. The dust tail 
hypothesis can probably be ruled out from geometrical 
considerations, but in any case it should have been 
observed by the more sensitive rocket instruments. For 
a neutral constituent to form the tail, it would have 
to be strongly influenced by radiation pressure yet have 
a very long lifetime, which is highly improbable at such 
a small sun-comet distance. The remaining candidate is 
an ion, with C"^ as the obvious choice, since it has a 
resonance multiplet in the bandpass (1335 A) and the sun 
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Vo = I km sec-' 




Fig. 6. Intensity contours for 01 1304 A computed for 
heliocentric velocities of 0 and -22.5 km 
sec"l. The outflow velocity of oxygen atoms 
is 1 km sec"l in (a) and 2 km sec"i in (b). 
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emits the mult ip let strongly. The fact that no ion 
tail was observed from the rockets after perihelion is 
explained by considering the motion of an ion created 
in a moving plasma. In the rest frame of a (perfectly 
conducting) plasma there is only a magnetic field, so 
the ion velocity consists of its original component 
parallel to the field plus a circular component due to 
its orbit around the field line. In the frame of the 
sun, the radial velocity of the ion will then consist 
of the sum of the radial velocity of the comet and a 
cycloidal component varying between zero and twice the 
radial plasma velocity (which can be as much as lO^ km 
sec"^ outward if the ion is in the solar wind). Thus 
before perihelion the ions spend considerable time with 
a small radial velocity and can resonantly scatter the 
solar line. After perihelion the radial velocity of the 
ions can never be less than that of the comet, so at 
the time of the rocket observations, when the ions 
were always Doppler shifted completely out of the solar 
lines, the ionized carbon was unobservable. 


CO Fourth Positive System 

Fluorescence of CO in the fourth positive system 
(A^tt . xls+) can be excited by solar radiation shortward 
of 1544 A. The g factors for the strongest bands are 
shown in Fig. 7, using oscillator strengths and 
branching ratios given by Mumma et al. (1971) and the 
continuum solar flux of Rottman Xprivate communication; 
see Donnelly and Pope 1973). Since excitation is from 
the continuum there is no variation in the g-factors 
with r. The g- factor for the total fourth positive 
emission between 1440 and 1670 A is simply the sum of 
all of the values shown in Fig. 7, 1.24 x 10"° photons 
sec"l- mol"^. The flux at the earth can be written as 

F 

4tt A'^ 


where Q is the production rate in molecules sec and 
the product g"^ is independent of heliocentric distance 
and may conveniently be evaluated at 1 A.U. In Table II 
we give the values of g and t for C, 0 and CO at the 
time of the rocket observations of Comet Kohoutek. If 
CO were in fact the parent of the observed carbon we 
would expect " Qq ~ t ^CO^ since the CO photodissociation 
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Fig. 7. g-factors for the strongest bands of the CO 
fourth positive system. 
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Table II 


Excitation Parameters for r = +55 km sec 


Relative 

Species g t gx Signal 

- 1-1 1 
(photons sec" mol’ ) (sec) (photons mol’ ) 


Cl X 1657 A 

2. 5 X 10"® 

2. 5 X 10^ 

6. 25 

.49 

OI X 1304 A 

3. 4 X lO’"^ 

4. 0 X 10® 

1. 36 

.28 

CO (1440-1670 A) 

1. 2 X 10'^ 

6. 9 X 10® 

0. 86 

.23 


(McElroy and McConnell 1971) and photoionization 
(Siscoe and Mukherjee 1972) rates are nearly equal. The 
table also gives the relative signal to be expected in 
this case, taking into account the response of a Csl 
photocathode, and it is seen that the three constituents 
give comparable contributions to the signal. Of course 
the primary source of oxygen, the photodissociation of 
OH, will give a larger 1304 A signal. 


It should be noted that detection of the strongest 
CO bands was just beyond the sensitivity of both rocket 
instruments, but that an order of magnitude improvement 
in either the instrument capability of the brightness 
of the comet at the time of observation would have 
enabled us to determine whether or not CO is the parent 
molecule. The values of derived from the data of 
Feldman £t al . (1974) have been revised using the g-factor 
discussed above, and are given in Table III. The upper 
limit on Q^q is about a factor of 3 higher than would be 
expected from the derived Q^. 


The (9, 0) fourth positive band at 1301 A is nearly 
coincident with the 01 line at 1302.17 A, raising the 
possibility of strong fluorescence in the v' = 9 
progression. However, the absorption oscillator strength 
is smaller than that for the (1, 0) band by roughly the 
same amount that the solar fliix in the line is greater 
than the continuum at 1510 A so that the absorption g- 
factors are of comparable magnitude. The rotational 
development of the band extends several A so that the 
overlap of the band with the relatively narrow 01 line 
is less than 10%, as shown in Fig. 8 for T = 100®K. 

The maximum g-factor obtainable for any one band in this 
progression is < 1.0 x lO"^ which is considerably smaller 
than the values given in Fig. 7. A plot of the^g-factor 
for the (9, 2) band at 1378 A as a function of r is shown 
in Fig. 9. 


Ho Lyman Band System 

Solar HI LyP can also produce fluorescence in the 
v' = 6 progression of the Lyman band system of Ho with 
the strongest emission in the (6, 13) band at 1608 A 
(Feldman and Fastie 1973). The PI line lies at 1025.935 A 
while the wavelength of the LyP line center is at 
1025.72 A, s§ that once again the g-factor is a sensitive 
function of r. The dependence on r is the same as the 
dashed curve in Fig. 4, except that the velocity axis 
is shifted so that the peak g-factors occur at values of 
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Table III 


SUMMARY OF RESULTS FROM AEROBEE 26. 023 


Species 

_F 

g 

T 



(photons sec ^ cm 

(sec mol 

(sec) 

(sec S 

OI 1304 A 

120 + 40 

3. 4 X 10'^ 

4. 0 X 10® 

2, 1 X 10^^ 

Cl 1657 A 

140 ± 50 

2. 5 X 10'^ 

2. 5 X 10® 

0. 5 X 10^^ 

OH 3090 A (0, 0) 

3100 + 100 

1. 2 X 10"^ 

1. 0 X 10® 

0. 6 X 10^^ 

CO 1510 A (1, 0) 

< 15 

2. 2 X 10"^ 

6. 9 X 10® 

< 2, 7 X 10^® 


RELATIVE INTENSITY 



Fig. 8. 


/ ■ 


Line strengths in the (9, 0) band of the CO 
fourth positive system showing the overlap 
with the solar 01 1203 A line. 
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r of +15 and +110 km sec . In this case, significant 
observations can be made only in the post-perihelion 
period. The maximum g-factor for the (6, 13) band is 
temperature dependent and has a value of 1.6 x 10“ 7 
sec"l mol“l for T = 300*K at 1 A.U. Using this value, 
the data of Feldman ^ al. (1974) give an upper limit 
to the H 2 production rate an order of magnitude smaller 
than the production rate of atomic hydrogen. 


SI A1300 A 


Finally, we consider the possibility of an 
accidental coincidence between a strong solar emission 
line and the resonance transition of a minor species. 

For example the sulfur multiplet near 1300 A contains a 
- ^Po line at 1302.34 which can resonantly scatter 
the 01 1302.17 A line over a wide range of heliocentric 
velocities. For this particular sulfur line the cooling 
time of the J=1 level is short (« 103 sec), so that 
resonance scattering from it is highly unlikely. However 
in the general case of wavelength coincidence, the 
identity of the scattering species would have to be 
determined either by means of high resolution spectros- 
copy or from a measurement of intensity over a large 
range of values of r. No unidentifiable features were 
found with either of the rocket spectrometers. 


Future Observations 

At present, the major uncertainty in the interpreta- 
tion of cometary ultraviolet emission features lies in 
the value of the species lifetime. Many authors have 
evaluated photoionization and photodissociation lifetimes 
but these are based largely on theoretical cross-sections 
In addition, the effect of solar wind ionization is not 
understood since the extent of solar wind penetration 
into the coma is unknown. Photographs of the cometary 
coma in monochromatic ultraviolet radiation, similar to 
the Lya photographs of Opal et al. (1974), if obtained 
at high enough angular resolution, should provide 
reliable values of the scale length and consequently of 
these lifetimes. Taken together with ultraviolet spectro 
photometric observations of the species discussed above 
over a large range of heliocentric velocities and 
distances (presumably from an earth orbiting satellite), 
it should be possible to determine accurately both the 
total gas production rate in the comet and its dependence 
on heliocentric distance. 
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DISCUSSION 


Voice : I just want to know if the values for production rates supersede the 
ones published in Science. 

P. D. Feldman : Yes. They will appear in the proceedings of this confer- 
ence. I don't think they are going to change. At least, we are not going to do 
any more work on it. 

I think what is called for now are some more observations. 


